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This  report  presents  the  results  obtained  irom  applying  an 
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buried  in  a stronger  off-azimuth  interfering  event.  Four  mixed-event  simu- 
lations with  various  adapt!  ve -filter  operational  specifications  are  performed  in 

this  report. 
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SECTION  I 
INTRODUCTION 


A.  PURPOSE  OF  TIRS  STUDY 

The  first  objective  of  the  signal  estimation  task  in  the  VELA 
Network  Evaluation  and  Automatic  Processing  Research  program  (Project 
VT/5705)  is  to  investigate  -md  report  on  new  techniques  for  optimally  detecting 
the  arrival  of  a second  seismic  event  buried  in  the  coda  of  a first-arriving 
event  at  various  signal-to-noisc  ratios. 

This  report  deals  with  results  obtained  from  operating  a maxi- 
mum likelihood  adaptive  beamforming  system  on  Korean  short-period  data  for 
mixed-event  separation. 


B.  DESCRIPTION  OF  KOREAN  SHORT-PERIOD  ARRAY 

Figure  1-1  depicts  the  geometrical  configuration  of  the  Korean 
Short-Period  Array,  which  is  a 19-element  hexagonal  array.  Table  1-1  gives 
the  array  sensor  locations. 

C.  DESCRIPTION  OF  THE  MAXIMUM  LIKELIHOOD  ADAPTIVE  FILTER 
ALGORITHM 

The  adaptive  filter  output  y(t)  at  time  t is  formed  by  applying 
a convolution  filter  to  each  channel  and  summing  the  outputs  of  all  filters: 

M N 


y<t)  = £ £ ai(j)  xi  (t-j) 


i = l j=-N 
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FIGURE  1-1 

KOREAN  SHORT-PERIOD  ARRAY  GEOMETRY 
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TABLE  1-1 

KOREAN  SHORT-PERIOD  SENSOR  LOCATIONS 

( Reference  Location  (37  27 ' 14"  N,  127  - 

15*  24"  E)  ) 

Sensor 

East  (X) 

North  (Y) 

Number 

Kilometers 

Kilometers 

-?..  34 
-1.  55 
-0.  26 
-0.05 
-2.  53 
-4.  32 
-4.  24 
-2.  85 

0.  13 
0.  86 
2.  90 

1.  93 
0.  81 


2.  24 
4.  35 
2.  37 


1. 61 

3.  48 
6.  52 
6.  00 

4.  34 
2.  60 
0.  85 

-1.  25 
-2.  57 
-1.  57 


1. 29 
3.  33 
5.  11 


where  a.IJ)  is  the  filter  weight  for  the  i-th  channel  at  a lag  of  j sample 
points,  x.  (t-j)  is  the  value  of  the  channel  i at  lime  t-j,  M is  the  number  of 

channels, ' and  2N+ 1 is  the  total  length  of  the  filter  in  points.  Prior  ,o 

forming  the  filter  output,  each  channel  is  time-shifted  to  time-align  energy 
arriving  from  the  desired  steer  direction. 

The  adaptive  filter  weights  are  updated  by  the  following 


algorithm  • 


new 


old 


a.(j)  - »,(1)  ♦ 2W 

i--  1 j = -N 


(1-2) 


where 


M 


x(t-j) 


= - x.(t-j) 

M Lu  i 


i=  l 


(1-3) 


and  M controls  the  adaptation  rate.  This  upda.e  algorithm  incorporates  the  maxi- 
mum  likelihood  constraints* 

,r  form,  the  adaptive-bcamforming  filter  update  equation 


In  vector 


may  be  written 


2 H XTA(X-X) 
yy(t+  At)  - A (t)  - t- 

X X 


(1-4) 


where  the  ouperscr.pt  T denotes  vector  transpose,  and  where  the  filter 
weight  vector  A,  the  data  vector  X.  and  the  beams.eer  output  vector  X are. 

respectively , 


1-4 


The  objective  of  maximum  likelihood  adaptive  beamforming 
is  to  reduce  the  average  squared  filter  output 

y2(t)  = (atx)  (xta)  = atxxta  {i“6) 

subject  to  a set  of  signal -preservation  constraints  on  the  filter  vec  tor  A. 
After  preshifting  the  input  channels  to  time -align  energy  from  the  look 
direction,  these  constraints  can  be  written 


M 

V' 


= -N 


-1,0,  1,  ...  N),  (1-7) 


where  the  constraints  d(j)  specify  a convolution  filter  having  the  desired  fre- 
quency response  on  a signal  from  the  steer  direction.  For  the  adaptive  beam- 
forming employed  in  this  study,  a fi?t  frequency  response  is  specified  by  set- 
ting 


d( j)  — 6jo  (j—  -N,  ...  ,-1,0,1,  ...  N) , 

where  8..  is  the  Kronecker  delta  operator 
i) 

1 if  i = j 

6..  = 

1J  0 if  i j . 

D.  ORGANIZATION  OF  REPORT 

Section  II  presents  the  results  of  mixed-event  simulations  which 
examine  the  possibility  that  adaptive  beamforming  can  detect  signals  buried  in 
off-azimuth  interfering  events  when  the  conventional  beamsteer  processor  can- 
not. The  simulation  process  is  performed  by  scaling  two  individually  recorded 
data  samples  and  summing  them  to  create  a composite  sample  used  in  the 
adaptive- filter  update  procedure.  The  adaptive  filter  is  then  applied  to  the 
two  data  samples  separately  to  form  two  separate  beam  outputs  which  are  sub- 
sequen'ly  added  to  create  a composite  beam  output.  The  detection  capability 
of  adaptive  processing  is  determined  from  a measurement  oi  the  amplitude  rise 
on  the  composite  trace  after  the  on-azimuth  event  arrival.  Ir  this  report,  the 
results  from  four  mixed-event  simulations  with  various  adaptive-filter  oper- 
ational specifications  using  Korean  short-period  data  are  presented.  Long- 
period  results  using  ALPA  data  have  been  included  in  a previous  report 
(Barnard  and  O'Brien,  1974), 

Section  III  presents  the  conclusions  of  this  study.  References 
are  given  in  Section  IV. 


d-8) 


d-9) 
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INTERFERING-EVENT  RESULTS 


A.  DISCUSSION 

The  ability  of  adaptive  beamforming  to  detect  an  on-azimuth 
event  buried  in  stronger  off-azimuth  interfering  events  has  been  successfully 
demonstrated  using  Alaskan  Long-Period  Array  data  [Barnard  ar.u  O'Brien, 
1974).  In  this  long-period  study,  maximum  likelihood  time-domain  adaptive 
beamforming  in  one  simulation  achieved  a 6 dB  amplitude  rise  after  the 
arrival  of  an  on-azimuth  signal  in  the  presence  of  an  interfering  event  30  dB 
above  the  target  signal  at  the  single-sensor  level.  For  Rayleigh  waves,  on 
the  average,  adaptive  beamforming  provided  an  array  gain  of  12  dB  (0.6  mag- 
nitude unit'-)  over  beamsteering.  In  addition  to  the  detection  improvement  for 
the  on-azimuth  signal  buried  in  the  off-azimuth  interfering  events,  adaptive 
beamforming  also  provided  a reasonably  accurate  estimate  of  the  on-azimuth 
signal's  maximum  peak-to-peak  amplitude  in  the  long-period  results,  where 
the  error  in  estimating  the  on-azimuth  signal  strength  was  less  than  0.  1 mag- 
nitude units.  These  promising  results  suggest  the  use  of  this  particular  algo- 
rithm for  P-wave  interfering-event  separation. 

The  next  subsection  describes  the  mixed-event  simulation  pro- 
cedure. Subsection  C discusses  t.>e  results  obtained  from  the  Korean  short- 
period  array  data.  These  results  are  limited  to  the  four  mixed-event  simula- 
tions conducted  during  this  contract  period  with  various  adaptive-filter  opera- 
tional specifications.  A brief  summary  of  the  results  is  presented  in  Subsection 
D. 


\ 


II- 1 
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B.  SIMULATION  PROCEDURE 

In  the  situations  where  two  wavetrains  from  two  separate  events 
pass  an  array  simultaneously,  these  wavetrains  will  overlap  to  form  mixed 
e\  ents.  Array  processing  techniques  attempt  to  separate  the  mixed  events  by 
steering  the  array  toward  the  target  signal  to  attenuate  the  off-azimuth  inter- 
fering-event energy.  However,  as  the  data  are  actually  recorded,  there  is  no 
evaluation  method  which  can  reliably  isolate  the  energy  contribution  from  eat  i 
of  the  two  mixed  events.  Hence,  the  most  realistic  method  to  evaluate  array- 
processing  techniques  such  as  beamsteering  or  adaptive  bcamforming  is  to 
rely  on  a simulation  procedure  in  which  two  data  samples,  each  containing 
one  signal,  are  added  to  create  a composite  sample  including  the  overlapped 
wavetrains  of  the  on-azimuth  signal  and  the  off-azimuth  interfering  event. 

Such  a simulation  can  be  designed  t > include  a number  of  parameters  through 
which  various  physical  situations  can  be  studied.  In  specifying  the  parameters 
for  a mixed-event  simulation,  the  relative  signal  strength  of  the  two  events, 
the  incoming  directions  of  the  two  signals,  and  the  time  separation  of  their 
arrivals  can  be  controlled  to  illustrate  the  real-world  situatiors. 

Figure  II- 1 presents  a schematic  diagram  for  the  simulation 
procedure.  As  shown  in  the  figure,  the  basic  mixed -event  simulation  procedure 
rums  two  data  samples:  one  contains  the  on-azimuth  signal  and  the  other  in- 

cludes the  off-azimuth  interfering  event.  Before  summing  the  two  samples  to 
form  a composite  sample,  each  input  channel  of  the  two  samples  passes 
through  the  same  zero-phase  prefilter.  The  on-azimuth  signal  sample  is 
time-shifted  to  time-nlign  the  target  event.  The  interfering  event,  on  the  other 
hand,  can  be  shifted  from  its  true  azimuth  and  apparent  P-wave  velocity  to 
the  desired  incoming  azimuth  and  veloc;'y  before  bcamforming  so  that  the 
interfering  event  will  appear  to  come  from  a fictitious  location  in  the  simulation. 
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BLOCK  DIAGRAM  OF  INTERFERING-EVENT  SIMULATION  PROCEDURE 


The  shifting  of  the  interfering  event  azimut  \ and  ve’ocity  Is  accomplished  by 
re-adjusting  the  time  delays  for  the  input  channels  of  the  interfering-event 
sample.  Following  the  formation  of  the  composite  sample,  a time-domain 
adaptive  multi-channel  filter  is  designed  from  it.  Initially,  a beamsteer 
filter  is  used,  and,  after  initialization,  the  adaptive  filter  updates  its  filter 
weights  according  to  the  update  equation 


A(t+  At)  = A + 


2MXTA  (X-X) 
T 

X X 


The  computer  software  package  can  simultaneously  update  the  filter  sets  for 
as  many  as  twenty  adaptive  convergence  rates  fiin  a single  processing  run. 

The  adaptive  filter  sets  are  subsequently  applied  to  each  individual  data  sample 
to  form  separate  beam  outputs  for  the  on -azimuth  signal  and  the  off-azimuth 
interfering  event.  Finally,  the  ABF  beam  for  the  on-azimuth  signal  and  the 
ABF  beam  for  the  off-azimuth  interfering  event  are  combined  to  form  the 
composite  ABF  trace.  The  time-shift-  and-  sum  beams  for  the  on-azimuth 
signal  and  the  off-azimuth  interfering  event  are  formed  similarly  and  subse- 
quently are  added  to  form  the  beamsteer  composite  trace.  The  composite 
traces  in  the  simulation  are  the  only  beam  outputs  which  are  available  in 
real-world  situations.  The  evaluation  of  adaptive  processing  performance 
relative  to  beamsteering  is  based  on  a comparison  of  the  composite  traces  for 
the  time-shift-and-sum  and  adaptive  beam  outputs.  However,  in  order  to 
facilitate  an  understanding  of  the  beamforming  processes  involved  in  the 
mixed-event  si'uations,  the  individual  sample  beams  for  the  on-azimuth  signal 
and  the  interfering  event  are  plotted  along  with  the  composite  trace  for  the 
beamsteer  output  and  the  specified  adaptive -beam  outputs.  Evaluation  of 
adaptive  beamforming  performance  relative  to  beamsteering  is  accomplished 
on  the  basis  of  the  composite-trace  beams. 
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1.  C.cneral 

The  principal  goal  of  this  subsection  is  to  examine  the  detection 
performance  of  adaptive  beamforming  relative  to  beamsteering  for  short-period 
interfering  events  simulated  by  adding  two  suitably  preprocessed  data  samples 
from  the  Korean  Short-Period  Array.  A secondary  objective  is  to  evaluate  the 
reliability  of  the  peak-to-peak  amplitude  on  the  composite-sample  adaptive 
filter  output  as  an  estimate  of  the  on-azimuth  signal's  peak-to-peak  amplitude. 
Because  this  is  the  first  application  of  adaptive  multichannel  filtering  to  Korean 
Short-Period  Array  data,  this  evaluation  also  provides  an  opportunity  to  gain 
an  understanding  of  the  optimal  choice  of  processing  parameters  such  as  the 
prefilter  applied  to  the  data  before  beamforming  and  the  azimuthal  separation 
between  the  interfering  event  and  the  target  event.  In  most  cases,  the  adaptive- 
filter  length  is  31  points  per  channel.  A shorter  filter  length  has  merit  in  on- 
line system  operation  if  it  does  not  substantially  degrade  ABF  performance. 
Therefore,  other  shorter  adaptive  filter  lengths  such  as  a 15-point  filter  length 
and  a 7-point  filter  length  are  also  tested.  At  various  points  in  this  subsection, 
four  events  are  used  to  simulate  mixed  events.  Table  II  — 1 lists  the  pertinent 
data  for  these  events.  Each  of  these  events  is  considerably  stronger  than  the 
background  noise  level.  Signal  similarity  varies  from  event  to  event  and 
from  passband  to  passband.  It  will  be  discussed  as  the  need  arises.  The  con- 
vergence rate  M for  the  adaptive  filter  upda-e  is  allowed  to  vary  in  each  pro- 
cessing simulation  in  order  to  observe  the  effect  of  the  convergence  rate  on 
adaptive  processing  detection  yain  relative  to  beamsteering. 

2.  Prefilters  Used 

Since  both  the  signal  and  noise  spectra  for  the  Korean  Short- 
Period  Array  data  are  still  not  well  understood  at  this  time,  a number  of 
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prefiltcrs  have  been  designed  for  the  mixed-event  simulations,  figures  II-2 
through  II - 6 show  prefilter  responses  in  dB  for  the  five  prefilters  employed 
in  this  study.  Figure  II - 2 displays  the  response  for  the  prefilter  with  a 1.0- 
2.0  Hz  passband  centered  at  1.  5 Hz.  The  .esponse  at  0,0  Hz  is  below  -120  dB 
so  that  the  d.  c.  bias  can  be  removed.  Figure  II - 3 shows  the  response  of  a pre- 
filter  with  the  lower-frequency  passband  0.  5 - 1.5  Hz  centered  at  1.0  Hz. 

The  corresponding  response  foi  a prefiltcr  with  a wider  passband  (1.0-3.  0 Hz) 
is  illustrated  in  Figure  II-4.  Figures  II - 5 and  II-6  display  the  response  of  two 
prefilters  with  0.  5-1.0  Hz  and  1.  25-1.75  Hz  passbands,  respectively. 


3.  First  Mixed-Event  Simulation 

To  simulate  the  first  mixed  event,  the  data  samples  for  the  Banda 
Sea  event  and  for  the  event  from  off  the  coast  of  Hokkaido  arc  summed.  The 
formcrevent  is  from  day  119  and  serves  as  the  on-azimuth  signal,  while  the 
latter  event  is  from  day  169  and  serves  as  the  off-azimuth  interfering  event. 

Start  times  for  the  two  data  samples  were  chosen  so  that  the  day  110  on-azimuth 
event  arrives  about  30  seconds  later  than  the  day  169  interfering  event.  The 
interfering-event  azimuth  was  shifted  from  64.  30°  to  357.  16°  to  achieve  a 180 
azimuthal  separation  from  the  target  event,  which  has  a 1m.  16  arrival 
azimuth  relative  to  the  Korean  Array.  The  P-wave  velocity  for  the  day  169 
intci fering  event  was  also  shifted  from  the  9.08  km/sec  true  apparent  velocity 
to  the  13.86  km/ sec  target  event  velocity.  Twelve  sites  (all  but  sites  1,2,  5,  9, 
10,  11  and  15)  constitute  the  channel  inputs  to  the  beamformers.  The  event 
separation  level  is  adjusted  by  scaling  either  the  on-azimuth  signal  or  the 
off-azimuth  interfering  event 

In  this  first  mixed-event  simulation,  the  data  traces  from  the 
day  119  event  were  scaled  so  that  the  averaged  peak-to-peak  amplitude  of  the 
target  signal  at  the  single-sensor  level  is  12  dB  less  than  the  corresponding 
averaged  amplitude  of  the  day  169  interfering  event. 
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Figure  II-7  illustrates  the  beamsteer  output  traces  for  the 
first  mixed-event  simulation.  The  upper  trace  is  the  beamsteer  outrut  for 
the  day  119  data  sample  with  the  on-azimuth  signal.  The  middle  trace  is  th 
beamsteer  output  for  the  day  169  data  sample  with  the  interfering  even,  shifted 
to  a 357.  16°  azimuth.  The  bottom  trace  is  the  composite-sample  beamsteer 
output  formed  by  summing  the  upper  and  middle  traces.  The  input  traces  for 
the  day  119  even,  were  scaled  to  a level  12  dB  below  that  of  the  day  169  inter- 
fering even,  so  that  the  composite  trace  for  the  adaptive  bcamformer  has  a 
6 dB  amplitude  rise.  Figure  II - 8 shows  the  adaptive  beam  outputs  for  the 
convergence  rate  H = 0.  3.  Signal  similarity  for  the  day  169  interfering 
event  is  poor  and  the  relative  channel  strength  varies  with  time.  This  event 
lasts  considerably  longer  than  a typical  P-wave  event.  The  amplitude  rise 
on  the  ABF  composite  trace  is  6 d!3,  while  the  beamsteer  output  yielded 

3.  8 dB  for  the  corresponding  pcak-to-peak  measurement.  In 

. ■ t o ,iri  nvrr  boamstoering  on  the  basis  of  pcak-to- 

best  ABF  improvement  is  d.  d all  over  ueam&i  ^ s 

peak  amplitude  measurements  on  the  composite  traces. 

In  Figure  II-8,  the  composite  trace  has  a 6 dB  borderline  do- 
lection.  Once  the  detection  is  declared,  the  next  mosj  important  consideration 
m evaluating  the  adaptive  bcamforming  performance  is  the  accuracy  with  which 
the  pcak-to-peak  amplitude  on  the  composite  trace  reflects  the  peak-to-peak 
amplitude  of  the  on-azimuth  signal.  Measurement  of  the  peak-to-peak  ampli- 
tudes on  the  upper  trace  and  composite  trace  for  the  on-azimuth  signal  window 
yields  the  same  amplitude,  as  shown  in  Figure  II-8.  Therefore,  the  amplitude 
on  the  composite-trace  adaptive  beam  for  this  mixed  event  furnishes  an  accurate 
bodywave  magnitude  estimate  for  the  detected  on-azimutl,  signal  if  the  adaptive 
resulting  apparent  period  of  signal  remains  unc.ianged. 

As  mentioned  earlier,  adaptive  beamforming  yields  only  2.2  dll 
improvement  over  beamsteering  on  the  basis  of  peak-to-peak  amplitude  measure- 
ments on  the  composite  traces  in  Figures  11-7  and  11-8.  However,  it  does  not 
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necessarily  follow  that  the  beamsteer  output  will  be  able  to  produce  a 6 dB 
amplitude  rise  for  an  on-azimuth  event  with  a bodywave  magnitude  0.  11  units 
higher.  To  demonstrate  this  fact,  the  day  169  on-azimuth  event  was  re- 
scaled so  that  the  beamsteer  output  on  the  composite  trace  is  able  to  achieve 
the  same  6 dB  amplitude  rise  on  the  basis  of  peak-to-peak  amplitude  measure- 
ments over  the  same  time  windows.  For  this  reason,  Figure  II-9  displays 
the  beamsteer  output  traces.  The  composite  trace  in  this  figure  has  the  same 
6 dB  amplitude  rise  as  the  corresponding  adaptive  trace  in  Figure  II- 8.  But 
the  off-azimuth  event  input  channels  in  Figure  II - 9 are  only  7.  6 dB  above  the 
on-azimuth  event  at  the  single  sensor  level.  Hence,  comparing  the  single 
sensor  levels  for  Figure  II- 8 and  those  for  Figure  II- 9,  adaptive  beamforming 
actually  yields  4.  4 dB  improvement  over  beamsteering  in  this  mixed-event 
simulation.  Figure  II- 10  shows  the  corresponding  adaptive  beam  output  at  the 
7.6  dB  event-separation  level,  where  the  amplitude  rise  on  the  composite  trace 

is  8.  5 dB. 

At  the  7.6dB  event  separation  level  used  to  produce  the  results 
of  Figures  II- 9 and  II-  1 0 , the  composite -trace  adaptive  filter  output  again  gives 
a more  accurate  estimate  of  the  event  strength.  For  the  beamsteer  output  on 
the  composite  trace  in  Figure  II-9,  the  peak-to-peak  amplitude  over  the  on- 
azimuth  signal  window  is  2.  2 dB,  or  0.  11  magnitude  units,  higher  than  that 
of  the  upper  trace  due  to  the  energy  from  the  interfering  event.  In  the  case  of 
the  adaptive  filter  output,  the  peak-to-peak  amplitude  on  the  composite  trace  in 
Figure  II- 10  is  about  1.  1 dB  lower  than  the  corresponding  amplitude  on  the 

upper  trace. 

To  examine  adaptive  beamforming  performance  at  various 
convergence  rates  for  this  mixed-event  simulation,  Table  II-2  tabulates  the 
amplitude-rise  results  and  the  ratios  of  the  composite-trace  amplitude  to  the 
corresponding  upper-trace  amplitude  over  the  on-azimuth  signal  window  for 
the  various  beams.  The  results  in  this  table  are  from  the  output  beams  at  the 


I 
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12  dB  event  separation  level.  Some  of  these  results  have  been  discussed  pre- 
viously in  comparing  the  detection  thresholds  for  beamsteering  and  adaptive 
beamforming.  At  low  convergence  rates,  adaptive  beamforming  suppresses 
the  on-azimuth  signal  more  than  the  off-azimuth  interfering  event.  Therefore, 
the  amplitude  rise  values  for  the  adaptive  convergince  rates  less  than  0.  1 are 
slightly  lower  than  those  of  the  beamsteer  output,  and  the  beamsteer  amplitude 
rise  is  equal  to  that  of  adaptive  beamforming  at  the  convergence  rate  /i  = 0.05. 
For  the  convergence  rates  greater  than  0.  1,  adaptive  beamforming  is  able  to 
suppress  the  interfering  e<ent  more  than  the  on-azimuth  signal.  T lis  fact 
results  in  the  ABF  improvements  over  beamsteering.  The  maximum  amplitude 
rise  in  the  table  is  6 dB  at  the  convergence  rate  M = 0.3.  For  the  convergence 
rates  greater  than  0.  3,  the  interfering  event  amplitude  remains  about  the  same, 
but  the  on-azimuth  signal  amplitude  is  lower  than  that  of  the  adaptive  beam  for 
M = 0.3.  As  a result,  the  amplitude  rise  is  lower.  The  ratio  of  the  composite 
trace  signal  amplitude  to  the  upper-trace  amplitude  for  the  time-shift-and-sum 
beams  is  3.  <■  IB.  For  adaptive  beamforming,  thiu  ratio  is  greater  than  3.4  dB 
for  the  convergence  rates  less  than  or  equal  to  0.01  and  is  equal  to  zero  for  the 
adaptive  beams  at  M = v.  3,  where  the  amplitude  rise  is  the  highest.  At  the 
adaptive  convergence  rates  greater  than  0.  3,  the  composite-trace  signal  is 
attenuated. 

Various  ABF  operational  specifications  were  tested  for  this 
mixed-event  simulation  using  the  events  from  days  119  and  169.  These  speci- 
fications include  the  prefilter  passbands  and  the  number  of  adaptive  filter 
weights  per  channel.  Table  II - 3 tabulates  the  amplitude  rise  on  the  composite- 
trace  for  the  beamsteer  and  adaptive  beam  outputs.  The  on-azimuth  signal 
attenuation  as  well  as  the  interfering  event  suppression  by  adaptive  beamforming 
relative  to  those  of  beamsteering  were  also  included  in  the  next-to-the-last  two 
columns  of  the  table. 
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In  the  0.  5-1.  5 Hz  passband,  the  waveform  similarity  for  the 
on-azimuth  event  of  day  119  is  poor  and  the  relative  strength  of  the  channels 
varies  with  time.  Therefore,  both  the  beamstecr  and  adaptive  amplitude  rise 
values  arc  the  poorest  among  the  passbands  shown  in  the  table.  The  A BP 
improvement  relative  to  beamsteering  is  also  the  poorest  for  this  mixed-event 
study.  Two  things  occurred  with  this  passband.  One  is  that  adaptive  beam- 
forming was  not  able  to  suppress  the  interfering  event  at  any  of  the  convergence 
rates  specified.  The  other  is  that,  at  the  higher  convergence  rates,  mutual 
cancellation  is  responsible  for  the  negative  results  relative  to  beamsteering 
at  the  higher  convergence  rates. 

In  one  case,  a 7-point  adaptive  filter  length  was  used  for  the 
l.  25-1.75  Hz  passband.  At  the  one -decisecond  sampling  interval,  the  7-point 
filter  length  covers  one  cycle  at  the  center  frequency  for  that  passband.  In 
comparison  with  the  results  in  the  next  row  of  the  table,  tlm  0.6  dB  increase 
in  the  ABF  gain  is  probably  not  enough  to  justify  using  that  adaptive  filtei 
length.  However,  the  best  convergence  rate  for  this  filter  length  is  H 0.4 
in  contrast  to  M = 0.0  5 for  the  31 -point  filter  length,  which  had  mutual  can- 
cellation at  the  higher  convergence  rates.  If  operating  at  the  high  converger  « e 
rates  only,  this  is  worth  noting. 

In  concluding  this  first  mixed-event  simulation,  the  adaptive 
amplitude-rise  results  are,  in  general,  about  2.0  dB  better  than  for  beam- 
steering.  The  waveform  similarity  and  the  time  stability  of  the  relative 
channel  strength  are  also  poor  for  the  day  169  interfering  event.  The  fact 
that  the  envelope  chape  for  the  input  channels  varies  significantly  across  the 
array  makes  it  difficult  for  the  ABF  processor  to  suppress  the  interfering  event. 
In  terms  of  the  ABF  detection  improvement  relative  to  beamsteering,  adaptive 
beamforming  can  reduce  the  bodywave  detection  threshold  by  about  0.  2 magni- 
tude units  from  the  beamstecr  level  for  this  mixed-event  simulation. 
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4.  Second  Mixed-Event  Simulation 

The  second  mixed-event  simulation  employs  the  same  Banda 
Sea  event  used  in  the  first  mixed-event  study  as  the  on-azimuth  signal  an  l uses 
the  Alaskan  event  from  day  MS  as  the  off -azimuth  interfering  event.  Originally, 
the  interfering  event  has  a 47.64°  azimuth  and  a 15.  1 km/sec  P-wave  velocity. 

In  this  simulation,  the  true  azimuth  of  this  even*  was  shifted  to  357.  16  and 
87.  16°,  while  its  true  velocity  was  shifted  to  that  of  the  on-azimuth  signal. 

The  start  times  for  both  samples  were  adjusted  so  that  the  on-azimuth  signal 
arrives  about  18  seconds  later  than  the  interfering  event.  This  interfering  event 
lasts  about  30  seconds,  considerably  less  than  the  first  one  from  day  169.  Four- 
teen sites  were  used  for  the  input  channels  (all  but  sites  9,  10,  11,  M,  and  19). 

Initially,  the  interfering-event  azimuth  was  shifted  to  3c-7.  16 
and  the  on -azimuth  signal  was  scaled  so  that  the  interfeiing  < vent  is  about 
18  dB  above  the  target  signal  at  the  single-sensor  level.  No  6 dB  amplitude 
rise,  was  achieved  on  either  the  ben  nsteer  as  the  ABF  output.  Subsequently, 
the  single-sensor  input  level  for  the  on-azimuth  signal  was  raised  to  1 2 dB  less 
than  that  of  the  interfering  event.  Using  a filter  length  of  31  points  per  channel 
at  this  event-separation  level,  adaptive  beamforming  was  still  unable  to  achieve 
a 6 dB  amplitude  rise.  However,  using  a 7-point-long  adaptive  filter  and 
shifting  the  interfering  event  azimuth  to  87.  16  or  90  separation  bc-t  ecu  tl  e 
two  events  .adaptive  beamforming  is  able  to  suppress  the  interfering  event 
enough  to  achieve  a 6 dB  amplitude  rise  at  tlu  event-separation  level  where 
the  interfering  event  is  16  dI3  above  the  on-azimuth  signal  at  tin  single-sensor 
level.  Figure  II-ll  displays  the  beamstcer  output,  where  the  amplitude  risi 
is  -1.6  dB  on  the  composite  trace.  The  corresponding  adaptive  beams  fo, 

U = 0.  5 are  shown  in  Figure  11-12.  A 6 dB  amplitude  lise  is  achieved  on  tlu 
composite-trace  beam.  At  this  high  convergence  rate,  adaptive  beamforming 
suppresses  the  interfering  event  11.5  dB  more  than  beamsteering,  but  degrades 
the  on-azimuth  signal  3.8  dB  relative  to  beamsteering.  This  fait  results  in 
the  7.7  dB  detection  gain  for  this  simulation.  The  signal  amplitude  on  tlu  AB1 
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composite  trace  is  3.  5 dB  lower  than  the  upper-trace  amplitude.  This  reduc- 
tion is  due  to  relatively  severe  mutual  cancellation  of  the  two  interfering  events. 

As  a matter  of  fact,  the  beamsteer  output  on  the  composite  trace  for  the  on- 
a/imuth  signal  was  2.  9 dB  higher  than  on  the  upper  trace  because  of  the  inter- 
fering-event energy.  Consequently,  the  actual  composite- trace  signal  degrad- 
ation by  the  ABF  processor  relative  to  beamsteering  is  about  1 dB.  In  order  to 
find  the.  detection  threshold  reduction  for  this  case,  the  event-separation  level 
was  re-adjusted  so  as  to  achieve  a 6 dB  amplitude  rise  on  the  beamsteer  output. 
Figure  II - 1 she  ws  the  time-shift-and-sum  beams  when  the  interfering  event  is 
5 dB  above  the  on-azimuth  signal  at  the  single-sensor  level.  In  other  words, 
the  threshold  for  a 6 dB  amplitude  rise  is  11  dB,  or  .55  magnitude  units  lower 
for  adaptive  beamforming  than  for  bcamsteering.  Figure  II- 14  shows  the  corres- 
ponding adaptive  beams,  where  an  18.8  dB  amplitude  rise  is  achieved  in  the  com- 
posite trace. 

Table  II - 4 tabulates  the  amplitude -rise  results  for  various  ABF 
operational  specifications  in  this  mixed-event  simulation.  Using  the  1. 0-2.0 
Hz  pass  band,  the  31 -[joint  adaptive  filter  length  at  180  azimuthal  separation 
\ i elded  a 2.  1 dB  improvement  relative  to  beamsteering.  A slightly  lower  im- 
provement value  is  obtained  using  a 7-point  filter  length.  Greater  improvement 
values  are  obtained  at  a 90°  azimuthal  separation  for  both  the  31  -point  and  7- 
point  filter  lengths  Using  the  narrower  1.  25-1.  75  Hz  passband,  the  detection 
gain  is  comparable  with  that  of  the  1.0-2.  0 Hz  passband  at  the  12  dB  event 
separation  level,  but  is  1.9  dB  less  for  the  larger  18  dB  event  separation  level, 
as  indicated  in  the  table.  A 6 dB  ABF  detection  gain  is  achieved  for  the  1.0- 
3.0  Hz  passband.  This  increased  ABF  improvement  can  be  attributed  to  mutual 
cancellation  on  the  time-shift-and-sum  beam.  This  case  is  an  exception.  For 
the  0.  5-1.  5 Hz  passband,  the  results  are,  in  general,  the  worst  among  the  pass- 
bands  shown.  A partial  explanation  for  these  results  is  the  poor  waveform  simi- 
larity and  the  time-varying  relative-  channel  strengths  for  the  on-azimuth  signal 
in  this  passband  (as  discussed  in  the  last  mixed-event  simulation.)  A 15-point 
adaptive  filter  length  was  also  specified  in  order  to  compare  its  adaptive  filter  per- 
formance with  that  of  the  31-point  filter  length.  In  this  case,  a 3,  5 dl  lower 
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Signal,  7 -Point  Adaptive  Length,  and  1. 0-2.0  Hz  Passband) 


detection  gain  relative  to  the  31 -point-filter  length  was  obtained.  1 he  180  a/.i- 

muthal  separation  for  this  mixed-event  from  days  119  and  1-15  yielded  only 
slightly  better  detection  gains  than  the;  previous  simulation  using  the  events 
from  days  119  and  1 69 , even  though  the  day  145  interfering  event  lias  better 
waveform  similarity  and  less  instability  in  relative  channel  strength.  However, 
at  a 90°  azimuthal  separation  between  the  intertering  event  and  the  on-azimuth 
signal,  an  adaptive-beamforming  amplitude  rise  11  dI3  (0.  55  magnitude  units) 
higher  than  that  of  beamsteering  was  achieved  with  the  same  two  events  using 
a 7-point  adaptive  filter  length.  Comparison  w'ith  the  0.  1 m^  unit  reduction 
achieved  with  a 180°  azimuthal  separation  and  a 31 -point  filter  length  suggests 
that  adaptive  beamforming  performance  varies  considerably  and  depends 
strongly  on  the  specific  operational  situation  encountered. 

9.  Third  Mixed- Event  Simulation 

The  third  mixed -event  simulation  combines  the  day  145  Alaskan 
event,  which  was  used  in  the  second  mixed  event  simulation  as  an  off-azimuth 
interfering  event,  and  the  day  125  event  from  the  Solomon  Islands  region.  The 
Alaskan  event  in  this  simulation  serves  as  the  on-azimuth  signal,  while  the 
Soloman  Islands  event,  with  a 143.69°  azimuth  and  a 52.85°  epicentral  delta 
relative  to  the  Korean  Short-Period  Array,  serves  as  the  off-azimuth  inter- 
fering event.  The  interfering  event  azimuth  was  shifted  to  227.64  and,  in 
some  cases  to  1 37.64°  to  achieve  180°  and  90°  azimuthal  separation  between 
the  two  events,  while  its  apparent  P-wave  velocity  was  shifted  to  15.  10  km/ sec, 
which  is  the  on-azimuth  event  P-wave  velocity.  The  start  times  for  the  two 
events  were  adjusted  so  that  the  on-azimuth  signal  arrives  20  seconds  later 
than  the  interfering  event.  Two  array  configurations  were  employed:  the 

7-element  inner-ring  array  and  the  19-element  field  array. 

Using  the  7 -site  inner  ring  array  (sites  1 through  7)  and  with 
the  interfering-event  azimuth  shifted  to  227.64  , the  on-azimuth  signal  was 
initially  scaled  so  that  the  interfering  event  is  18.  8 dB  above  the  on-azimuth 
signal  at  the  single- sensor  level.  Neither  the  beamsteer  nor  the  adaptive 
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processor  aclncvcu  a o elli  anipl nude- rise  value.  In  orde  r to  compare  tin 
detection  thresholds  between  the  ABF  and  beamsteer  processors,  as  was  done 
previously  in  the  last  two  mixed-event  simulations,  the  on-azimuth  signal  was 
scaled  to  3.  1 dB  less  than  the'  interfi  ling  event  at  the  single-sensor  level  so  that 
a 6 dB  amplitude-rise  value  for  the-  ABF  beam  eoulei  be  achieved  Figure  11-19 
displays  the  beamsteer  output,  where  a 0.9  dB  amplitude  rise-  oe  curreel  on  the 
composite-trace  time-shift  nd-sum  beam.  The  cor  re-sponeling  adaptive  beams 
for  the  convergence-  rate  0.9  are  shown  in  Figure  1 1 - 1 6 , where  a mlB 

amplitude  rise  occurred  on  the  composite-trace  output.  Among  the-  tinm-shift- 
and-sum  beams  in  Figure  11-19,  the  composite-trace  on-azimuth  signal  is  1.9 
el  B 1 e'ss  than  on  the  upper  trace,  while  in  Figure  II- 16  the  coinposi  te- trace  on- 
azimuth  signal  amplitude1  is  0 .2  dB  less  than  that  on  the-  upper  trace.  Therefore, 
the  ABF  processing  gain  is  mostly  due-  to  interfering  e vent  suppression.  In 
order  to  achieve'  a 6 dB  amplitude  rise-  for  beamsteering,  the'  on-azimuth  signal 
has  to  he-  re-scaled  so  that  it  is  2.  9 dB  above-  the'  inte-rfe  ring  e-vent  at  the  single- 
sensor  1 eve-1.  Figure  II- 17  shows  the  beamsteer  output  for  this  case.  Because1 
the  input  channels  for  the-  on-azimuth  signal  are1  higher  than  those  evf  the-  uite-r- 
fering  event,  the1  composite-trace  on-azimuth  signal  is  essentially  equal  to  that 
on  the  uppe r trace.  Figure  11-18  displays  the1  corresponeling  aelaptive  beams, 
where  the1  amplitude  rise  on  the  composite-  traee*  is  17.  1 dB,  In  comparisem 
with  the  be-amstewr  detection  threshold,  aelaptive  beamforming  is  able  to  reduce 
the  detection  lcvc-1  b^  9.9  e!B,  e-epiiv  ulent  to  0.28  magnitude  units  for  this  rnixeel 
event  with  the  inner -ring  array. 

A number  of  case's  were1  processeel  for  this  mixeel -e-ve-nt  simula- 
tion. Table-  II  - 5 tabulates  the-  ampl  it  ude  - rise  re-suits'  using  the  7-site-  inner- 
ring  array.  In  two  cases,  the  inte  rfe-ring  -e-vent  azimuth  was  shifte-el  from  199.69 
to  137  69  to  achie-ve'  a 90  azimuthal  separation  he  tween  the  two  e vents.  In 

two  cases,  an  1 1 -point  filter  length  was  eompareel  with  a 31-point  filter  length. 

o 

None-  eif  tlu-se-  ])rovicle-el  a highe  r ABF  gain  titan  a 31-point  filter  with  a 180 
azimuthal  separation,  as  inelicateel  in  the’  table-.  In  addition  to  the  0.  9-1  9 11/ 
passband,  0.  9-1.0  Hz  anel  l. 0-3.0  Hz  oassbands  we  re  also  used.  Among 


1 


1 

I 


II-  31 


I PQ 
In  o c 3 
PQ  S o c 

< « ~ ns 

Q O 


m r- 
o O' 


r-  ^ 

00  lO 


•PJ  -M  cn  « ' 

Sh  c w ,2  P Mi 

H)  O O *.  5 c 

V ^ «s  £ r 

<u  w S*  '3  pq  <u 

c §*«  £ 

H 10  to 


r-  oo  in  in 

• • • • 

m oo  r~- 


O (M  O 

N iri  ifi 


d ns  (tj  •£  ns  i? 

u*  £ <=  3 * it  -5 


<U  1) 
Oh  2 1) 


tM  OH  ITl  00  o vO  O 

• • • • * * * 

r-  if  M oo  vO  o v£> 


in 

i co 


0^0 
cr;  z 

W W a 
in  Z u 
< 2 Im 
U M « 
io  W W 

^ h u 


•*-'  f ^ D c. 

w H CX  -Ti  <y 

^ cQ  £ a tn 

w < < ~ £ 


I I CQ 

£ O ry  QJ  * — 

£ O c TO  <u 

^ *-*  C 3 in 


PQ  1,1  <; 


~ ns  PQ 

S >i  "C 

£ ns  'r- 

> Q,  C 

W 4J  o 

IO  '-2 


« 5 c g S 
cq  < o cc  & 
U IX 


' nJ  2 

NX  H C 

" « ns  O 

< d a 


M ~ <u  in 

In  v To  m c t! 

CQ  X .h  o g .5 

< ;H  H 2 o 

^ ^ U ^ 


I U l 

o;  a;  cn  tj  *?? 

m ~ w G £ 

Ql  I 

*m  CL  .O  w 


m 

• • 
O -« 


O'  O 'O 

« • ■ 

o ^ o 

I *“H 


oo  oo  oo  m 


«J*  SO 


oo  oo  oo  m m M 


D 0 0 0 0 0 0 

o o o o o o o 

00  O'  00  O'  00  00  00 


Only  one  convergence  rate  was  used 


these  passbands,  the  0.  5-1.  5 Hz  band  produces  the  best  results,  as 
shown  in  Table  11-5.  In  this  simulation,  adaptive  beamforming  achieved  a 
substantial  improvement  relative  to  beam  teering.  This  improvement  can  hi 
attributed  mostly  to  the  poor  beamstcer  performance  at  that  steer  velocity 
and  passband.  Figure  11-19  shows  the  time-shift-and-sum  beam  pattern  for 
the  47.6°  steer  azimuth  and  the  15.  1 km/sec  velocity  at  1.0  Hz  for  the  7-site 
inner- ring  array  and  the  19-site  full  array.  The  7 -site  pattern  is  almost  cir. 
cular  with  unity  response  at  the  look  direction  and  has  a -6  dll  value,  in  con- 
trast to  -20  dB  for  the  19-site  pattern  in  the  interfering-event  separation.  It 
reflects  the  low  beamsteer  array  gain,  which  ranges  from  5.6  dB  to  8.4  dB 
for  the  2.4  dB  to  -18.8  dB  event  separations,  as  shown  in  I able  11-5. 


This  mixed -event  simulation  was  repeated  using  the  19  sites  ol 
the  full  array  as  input  channels.  Figure  11-20  shows  the  time- shift -r-  nil  - sum 
beams  for  the  0.5  -1.5  Hz  passband.  By  increasing  the  number  of  sites  to 
enlarge  the  aperture  of  arr,v;  , a beamsteer  array  gain  of  nearly  20  dB  was 
obtained.  In  contrast  to  the  substantial  ABF  improvement  relative  to  beam- 
steering achieved  by  using  the  7-site  array,  acaptive  beamforming  using  the 
19-site  full  array  provided  no  improvement  with  respect  to  the  beamsteer 
results.  Figure  11-21  shows  the  best  ABF  beams  at  the  convergence  rate 
M - 0.1,  where  the  amplitude  rise  on  the  composite  trace  is  only  comparable 
with  the  beamsteer  output.  However,  one  may  be  able  to  see  a gap  between 
the  two  event  wavetrains  on  the  adaptive  composite -trace  beam  in  Figure  11-21. 


At  the  higher  convergence  rates,  adaptive  beamforming  attentnates  the  on 


azimuth  signal  rapidly  so  that  the  amplitude  rise  decreases.  The  adaptive 
beamformer,  in  order  to  minimize  the  power  output,  creates  a waveform 
180°  t of  phase  with  the  upper  trace  on  the  interfering-event  beam  output 
and,  as  a result  produces  mutual  cancellation  on  the  composite  trace. 
Figure  11-22  illustrates  these  situations  for  the  ABF  beams  at  the  convergi 
rate  M = 0.  a.  In  this  mixed-event  simulation,  mutual  cancellation  appeal 
in  other  passbands  at  high  convergence  rates.  Table  11-6  tabulates  the 
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FIGURE  11-10 
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5,  Signal  Azimuth  47.  6°,  Signal  Apparent  Velocity  15.  1 km/sec,  Interfering 
Azimuth  227.6°,  Interfering  Event  18.  8 dB  Above  the  On-Azimuth  Signal 
19  Sites  and  0.  5-1.  5 Hz  Passband) 


amplitude  -rise  results  for  three  passbands  where  only  adaptive  convergent  c 
rates  greater  than  0.  1 were  specified  in  the  simulation. 

6.  Fourth  Mixed-Event  Simulation 

The  ’ast  mixed-event  simulation  was  performed  using  the  day 
145  event  from  the  south  of  Alaska  as  the  on-azimuth  signal  and  the  day  169 
event  from  off  the  coast  of  Hokkaido  as  the  off-azimuth  interfering  event. 

Both  the  events  have  been  used  in  the  previous  mixed-event  simulations. 

In  this  mixed-event  simulation,  fourteen  sites  were  used  (all  but  sites  1,2,5, 

9 and  15)  The  start  times  for  the  two  events  were  adjusted  so  that  the  on- 
azimuth  signal  arrives  55  seconds  after  the  arrival  of  the  interfering  event. 

In  the  previous  simulations,  where  adaptive  beamforming  out  performed  beam- 
steering, the  best  adaptive  processing  results  are  usually  at  high  convergence 
rates.  However,  in  this  simulation,  mutual  cancellation  occurred  at  high 
convergence  rates.  As  a result,  the  best  adaptive  processing  results  are 
only  comparable  to  those  of  beamsteering.  Table  II-7  tabulates  the  amplitude 
rise  results  for  the  various  processing  parameters  used.  The  31 -point  adaptive 
filter  length  produces  the  best  results  among  me  ttu’.'C  filter  lengths  specified 
for  the  0 5- 1 . 5 Hz  pas sband.  For  the  1. 0-2.  0 Hz  and  1 0 -3.  0 Hz  passbands, 
the  adaptive  gains  relative  to  beamsteering  are  only  comparable  to  those  of 
the  0.  5-1.  5 Hz  pas  sband,  but  both  the  adaptive  and  beamstcer  array  gains 
with  respect  to  the  single-sensor  signal -to -noise  ratio  for  these  passbands  are 
about  9.0  dB  higher  than  those  of  the  0.  5-1.  5 Hz  passband;,  as  shown  in  the 
table.  For  the  1.0-2.  0 Hz  and  1. 0-3.0  Hz  passbands,  the  beamstcer  composite- 
trace  amplitude  increased,  because  of  the  interfering-event  energy,  by  about 
i , y dB  relative  to  the  upper  trace.  This  fact  suggests  that  the  low  beamstcer 
array  gain  for  the  0.  5-1.  5 Hz  passband  is  due  to  the  poor  time- shift -and- sum 
beam  pattern  for  that  passband. 
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SUMMARY 


In  the  first  mixed-event  simulation,  the  adaptive  amplitude 
rise  is  about  2.  0 dB  higher  than  for  beamsteering  in  most  of  the  passbands 
processed  except  for  the  0.  5-1.  5 Hz  passband,  in  which  the  on-azimuth  signal 
similarity  is  very  poor.  The  adaptive  processor  can  reduce  the  detection 
threshold  0.  22  bodywave  magnitude  units  on  the  basis  of  the  mixed  event  single- 
sensor amplitude  ratios  at  which  a 6 dB  amplitude  rise  occurs  on  the  beamsteer 
and  adaptive  beam  outputs.  Adaptive  beamforming  is  able  to  suppress  the  off- 
azimuth  interfering  event  prior  to  the  on-azimuth  signal  arrival,  but  it  seems 
less  effective  in  suppressing  the  off-azimuth  event  after  the  on-azimuth  signal 
passes  the  array,  particularly  in  the  case  of  wider  passbands. 

The  second  mixed-event  simulation  employs  the  same  on-azimuth 

signal  as  the  first  mixed-events  and  uses  an  off-azimuth  interfering  event  with 

better  signal  similarity  and  a shorter  duration  than  the  first  interfering  event. 

The  results  are  only  slightly  better  than  those  of  the  first  mixed-event  simula- 

o 

tion.  However,  for  the  cases  where  the  interfering  event  was  shifted  to  a 90 
azimuthal  separation  relative  to  the  on-azimuth  signal  and  where  a 7-point 
adaptive  filter  length  was  used,  a 0.  55  mb  detection  threshold  reduction  occurred 
with  adaptive  processing  on  the  basis  of  the  single -sensor  signal -to -noise  ratios 
at  which  a 6 dB  amplitude  rise  occurred  on  the  beamsteer  and  adaptive  outputs. 

In  the  third  mixed-event  simulation,  the  processing  was  performed 
using  the  7-site  inner-ring  array  and  the  19-site  full  array.  For  the  cases  with 
the  7-site  inner-ring  array,  the  beamsteer  array  gain  relative  to  a single  sensor 
ranges  from  3.6  dB  to  8.4  dB,  while  the  corresponding  ABF  array  gain  is  from 
12  dB  to  20.4  dB.  For  the  cases  where  beamsteering  and  adaptive  beamforming 
achieve  a 6 dB  amplitude  rise,  adaptive  beamforming  reduces  the  detection  thres- 
hold by  6.  5 dB,  or  0.  28  m,  units.  For  the  cases  with  the  19-site  full  array  in 
the  same  mixed-event  simulation,  the  best  ABF  amplitude  rise  is  only  close  to, 
but  not  better  than,  that  of  the  beamsteer  output. 
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The  last  mixed-event  simulation  resulted  in  situations  where 
mutual  cancellation  occurred  on  the  adaptive  beam  outputs  at  high  convergence 
rates  and,  as  a result,  the  adaptive  detection  performance  is  only  comparable 
to  that  of  beamsteering.  The  average  adaptive  detection  gain  relative  to  beam- 
steering is  about  0.6  dB  for  the  five  cases  tried  in  this  simulation. 

In  this  mixed -event  study,  array  processing  gain  is  the  differ- 
ence between  interfering-event  suppression  and  signal  degradation.  For  time- 
shift-and-sum  beamforming,  the  processing  gain  is  a constant  for  a given 
array  configuration,  steer  velocity,  passband,  and  look  direction,  and  is 
independent  of  the  event-separation  level  between  the  mixed  events  involved. 
However,  in  the  case  of  signal  degradation,  the  composite-trace  signal  ampli- 
tude was  affected  by  addition  or  subtraction  of  interfering-event  energy, 
which  could  produce  inaccurate  magnitude  estimates.  For  example,  in 
Table  II-2,  the  composite-trace  beamsteer  signal  amplitude  is  3.4  dB  higher 
than  that  on  the  upper  trace.  This  higher  value  is  due  to  the  addition  of  the 
interfering  event  to  the  signal.  For  the  optimum  adaptive-beamforming  pro- 
cessor in  this  case,  both  the  upper  trace  and  the  composite  trace  have  the  same 
amplitudes,  as  shown  in  Table  II- 2,  where  their  ratio  is  0 dB  at  M = 0.3.  An 
a consequence,  the  time-shift-and-sum  beamformer  registers  a higher  on- 
azimuth  signal  magnitude,  and,  in  contrast,  the  adaptive  beamformer  provides 
an  accurate  measurement. 

The  results  presented  in  this  section  arc  obtained  from  a maxi- 
mum likelihood  adaptive  algorithm,  where  the  design  goal  is  to  minimize  the 
filter  output  power  subject  to  unity-response  constraints  in  the  beam  look  direc- 
tion. The  adaptive  performance  here  is  not  as  good  as  its  performance  for  long 
period  ALPA  data.  Better  detection  might  be  possible  with  a Wiener  time- 
domain  adaptive  algorithm,  where  the  design  goal  is  to  minimize  the  mean 
square  difference  between  the  filter  output  and  the  on- azimuth  signal. 
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SECTION  III 


CONCLUSIONS 


In  the  mixed-event  simulations  presented  in  this  report,  two 
data  samples,  one  with  an  on-azimuth  signal  and  the  other  containing  an  inter- 
fering event,  are  summed  to  form  a composite  sample  from  which  the  adaptive 
filter  sets  are  designed  and  applied  to  the  two  samples  individually.  The  out- 
puts of  the  two  samples  are  added  to  form  the  composite  beam  used  in  simu- 
lating the  actual  situations  encountered  in  practice.  Amplitude-rise  measure- 
ments are  made  by  taking  the  ratio  of  the  composite -trace  maximum  peak-to- 
peak  amplitude  after  the  signal  arrival  to  the  corresponding  amplitude  before 
the  signal  arrival. 


The  results  vary  from  case  to  case.  The  following  points  sum- 


marize the  results: 


In  the  first  mixed-event  simulation, where  the  interfering  event 
has  poor  waveform  similarity  and  time-varying  relative  ampli- 
tudes across  the  array,  a 0.  2 m^  detection  threshold  reduction 
with  respect  to  the  beamsteer  level  was  obtained  by  adaptive 
beamforming.  The  second  simulation  uses  the  same  on-azimuth 
signal  as  in  the  first  and  employs  an  off-azimuth  interfering 
event  with  better  waveform  similarity.  The  results  are  slightly 
better  than  those  of  the  first  simulation.  However,  with  a 7- 
point  long  adaptive  filter  and  with  a 90°  azimuthal  separation  be- 
tween the  two  events,  adaptive  beamforming  fas  able  to  reduce 

the  detection  threshold  by  0.  5 m,  units  from  the  beamsteer  level. 
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Using  the  7-site  inner-ring  array  (sites  1 through  7),  a detec- 
tion threshold  reduction  of  about  0.  3 magnitude  units  relative 
to  beamsteering  was  achieved  in  the  third  mixed-event  simula- 
tion. In  this  case,  the  adaptive  gain  can  be  attributed  to  the 
poor  tim  e - shift- and  - s um  be?  m pattern.  With  the  19-site  full 
array  in  this  same  simulation,  the  adaptive-processing  results 
are  only  comparable  to  those  of  beamsteering  because  mutual 
cancellation  occurs  at  high  convergence  rates.  The  last  mixed- 
event  simulation  only  produces  comparable  results  between  the 
adaptive  and  beamstecr  processors  due  to  the  occurrence  of 
mutual  cancellation  at  high  convergence  rates. 

The  bermsteer  array  gain  relative  to  the  single-sensor  level  is 
more  affected  by  addition  or  subtraction  of  interfering-event 
energy  than  that  of  adaptive  beamforming  in  this  mixed -event 
study.  As  the  event-separation  level  increases,  the  array  gain 
tends  to  increase.  For  on-azimuth  signal  magnitude  estimates, 
the  adaptive-beamforming  processor  seems  to  produce  more 
accurate  results  than  beamsteering. 

Among  the  various  passbands  used,  the  1-Hz-widc  prefilter 
passband  yielded  the  best  results.  The  adaptive  gain  for  various 
passbands  is  comparable  as  long  as  signal  similarity  does  not 
vary  significantly  among  the  passbands.  After  the  arrival  of 
the  on-azimuth  signal,  adaptive  beamforming  tends  to  suppress 
the  off-azimuth  interfering  event  less  with  the  wider  passbands. 

In  these  short-period  P-wave  mixed-event  simulations,  both  ‘no 
beamsteer  and  adaptive  proc.  ssors  did  not  perform  as  well  as 
in  the  long-period  mixed-event  simulation  using  ALPA  data, 
particularly  in  the  case  of  the  adaptive  processor.  .Signal  simi- 
larity for  the  short-period  P-wavc  signals  from  the  spring  1)7-1 
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Korea  data  is  not  as  good  as  for  the  Rayleigh  waves  processed 
in  the  long-period  study.  Therefore,  adaptive  beamforming 
performs  much  better  for  the  long-period  Rayleigh  waves  than 
for  the  short-period  P-waves. 


In  the  spring  of  197  3,  the  Korean  short-period  arra  data 
appeared  to  have  some  digitizer  problems.  Whether  these  problems  have  signifi- 
cant effects  on  the  results  presented  in  this  repori  is  still  not  known.  The  digiti- 
zer problems  seem  to  be  corrected  in  the  November  1974  data.  On  the  basis 
of  a limited  number  of  samples,  the  data  quality  during  this  period  seems  to 
have  been  greatly  improved  and  seems  to  have  better  signal  similarity  across 
the  array.  Repeating  this  work  with  the  November  1974  data  might  produce 
significantly  better  detection  results. 
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